A search for non-virialized axionic dark matter 
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Cold dark matter in the Milky Way halo may have structure defined by flows with low velocity 
dispersion. The Axion Dark Matter experiment high resolution channel is especially sensitive to 
axions in such low velocity dispersion flows. Results from a combined power spectra analysis of 
the high resolution channel axion search are presented along with a discussion of the assumptions 
underlying such an analysis. We exclude KSVZ axion dark matter densities of p > 0.2 GeV/cm 3 
and DFSZ densities of p > 1.4 GeV/cm 3 over a mass range of m a = 3.3^teV to 3.69/ieV for models 
having velocity dispersions of A/3 < 3 x 10~ 6 . 



The Cold Dark Matter (CDM) component of the Milky 
Way halo may consist of low mass (/xeV-meV) pseu- 
doscalar particles called axions, which are a consequence 
of the Peccei-Quinn solution to the strong CP problem[l|- 
3]. Axions of this mass range would have extremely 
weak couplings to standard-model particles as well as 
to both axionic and non-axionic dark matter, rendering 
them effectively collisionless0-[l]- Models of the struc- 
ture of halos consisting of such particles predict discrete 
flows with low velocity dispersion and high density at 
special locations called caustics 0-11 • In a Sikivie type 
detector [13, [Hj], a low dispersion axion flow would ap- 
pear as a narrow peak in the power spectrum. The Ax- 
ion Dark Matter experiment (ADMX) high resolution 
(HR) data acquisition channel is sensitive to such sig- 
nals, even more so than for signals exp ected from halo 
models lacking low dispersion flows[12j|. Sensitivity to 
these signals is further improved by combining data from 
successive power spectra. This gain in the signal to noise 
ratio (SNR) relative to a single-scan analysis, such as in 



the previous ADMX HR analysis found in Ref. [12J, al- 
lows for sensitivity to lower densities, though care must 
be taken to account for the change in frequency of the ax- 
ion signal due to the Earth's motion relative to the dark 
matter. Because any measurement of the axion's kinetic 
energy will be subject to Doppler shifts, both spectral 
broadening and signal modulation are expected. In this 
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paper we discuss the constraints on a combined-spectra 
analysis imposed by Doppler shifts, and present the re- 
sults of the first ADMX combined-spectra HR analysis. 

ADMX converts axions to detectable microwave pho- 
tons via the inverse Primakoff effect within a tunable, 
high quality factor (Q ~ 50,000) microwave cavity im- 
mersed in a strong magnetic field. An axion with mass 
m a and velocity /3 would create a photon with energy 
E ss m a c 2 + \m a c 2 f3 2 . A velocity dispersion A/3 yields 
spectral broadening of ^j- = 2^jr « 2/3 A/3. Axion con- 
version is expected to produce power in the cavity of 
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where V is the volume of the cavity, Bq is the magnetic 
field strength, and p a is the axion density in the detector. 
The mode dependent form factor of the cavity is given 
by 
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where E is the electric field of the tuned mode and Bo is 
the external magnetic field. The highest form factor ob- 
tainable (C ss 0.69) corresponds to the TMoio mode[l3[ 
prompting its use over higher order modes. The quality 
factors of the cavity and the axions are Q < 10 5 and 
Qa ^ 10 6 , respectively. Finally the axion to two photon 
coupling constant (7,177 is given by 
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FIG. 1. Diagram of the ADMX cavity and receiver chain. The power read out from the cavity is sent through both cold 
and room temperature amplification stages, is mixed down from radio frequencies to audio frequencies via two local oscillators 
(L.O.), and is ultimately saved to disk. 
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where a is the fine structure constant, f a is the Peccei- 
Quinn symmetry breaking scale, and g 7 is a model 
dependent constant. In the Kim-Shifman-Vainshtcin- 
Zakharov (KSVZ) model g 1 = -0.970 [3, and in the 
Dine-Fischler-Srednicki-Zhitnitskii (DFSZ) model g 1 = 
0.3600. 

The ADMX detector consists of a 1-m tall, 0.5-m di- 
ameter, copper plated, stainless steel, right cylindrical 
cavity kept at 1.8 K and placed in a 7.6 T magnetic field. 
The TMoio mode of the cavity is swept through a range 
of frequencies, corresponding to a range of possible axion 
masses, by moving a metal rod, oriented parallel to the 
cylinder axis and extending along the full length of the 
cavity, from the wall to the center of the cavity space 0]. 
For the experimental parameters of ADMX, Eq. [1] pre- 
dicts an expected signal power of order 10 -22 W. The 
expected SNR is given by the Dicke radiometer equation 



where P is the expected signal power, P/v is the noise 
power, b is the signal bandwidth, t is the integration 
time, ks is Boltzmann's constant, and Tjv is the total 
noise temperature which is equal to the physical temper- 
ature plus the noise temperature of the electronics. 

ADMX has two data acquisition channels which use 
the receiver and amplifier chain shown in Fig. [TJ Power 
from the cavity is critically coupled to the receiver which 
sends the signal through a directional coupler followed 
by two cold (~ 1.8 K) amplification stages and then one 
room temperature amplification stage. The cold am- 
plification stages are a DC Superconducting QUantum 
Interference Device (SQUID) amplifier [20j followed by 



two balanced GaAs heterostructure field-effect transis- 
tor (HFET) amplifiers. There is ~10 dB gain from the 
SQUID and - 34 dB combined gain from the HFETsHH- 
|23[ . At room temperature, the signal is amplified an ad- 
ditional 35 dB before being down-converted to a center 
frequency of 10.7 MHz and sent through an eight-pole 
crystal band pass filter. The signal is down-converted a 
second time to a center frequency of 35 kHz whereupon 
the medium resolution (MR) and HR channels diverge. 

The MR channel takes 10,000 scans over 80 seconds, 
computes a fast Fourier transform (FFT) and power spec- 
trum of each scan and averages the 10,000 spectra into 
one 400-point power spectrum with a resolution of 125 
Hz. The HR channel takes data for 23.8 seconds with 
a sampling rate of 80 kHz resulting in a maximum reso- 
lution of 42 mHz. Three such spectra are taken during 
the integration time of the MR channel. The HR data 
are saved in two formats. First the raw time series data 
are written to disk with no signal processing or averaging 
having been applied. Analysis of these data will be the 
subject of a future publication. Second, the three scans 
taken during the MR integration are Fourier transformed, 
averaged, and saved as a single 42-mHz bin- width-power 
spectrum. After acquiring data at a given resonant fre- 
quency, /o, the tuning rods are moved, shifting /o ~ 2 
kHz, and the data acquisition process is repeated. Be- 
cause the tuning rod steps are much smaller than the 
30-kHz bandwidth of the crystal filter and the resonant 
width of the cavity TMoio mode, the frequency coverage 
for successive scans overlap significantly. For both chan- 
nels, the total exposure time for a given frequency bin is 
~25 minutes. 

The MR channel has recently produced limits for both 
virializcd and non-virializcd axions with A/3 ~ 10~ 3 and 
< 2 x 10" 4 , respectively^. While the MR and the HR 
channels both look for non-virialized axions, the differ- 
ence comes in their sensitivity to models with differing 
velocity dispersions. The MR channel is equally sensitive 
to all signals with dispersions of < 10~ 4 , as all such sig- 
nals would occupy a single 125 Hz bin. The HR channel 
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FIG. 2. Monte Carlo simulations of non-virialized axion sig- 
nals imposed on real HR data. The broader peak (velocity 
dispersion of 2 x 10~ 4 ) is shown as an example of the narrow- 
est signal that the MR channel can resolve. An excess power 
of corresponds to power in that bin equal to the rms noise 
power at that frequency. 



potentially has a much smaller bin width, and therefore 
continues to gain in the SNR for signals with low veloc- 
ity dispersion. Figure [2] shows the expected signals in the 
HR channel for two simulated axion flows of equal den- 
sity but differing velocity dispersion. While both signals 
would be be seen in the HR channel, the signal having 
lower velocity dispersion stands dramatically above the 
noise. Note that at its highest resolution, the HR channel 
would have non-Gaussian noise. For a detailed discussion 
of the possible origins and the noise characteristics of low 
dispersion axion signals see Duffy et al.fl^. 

A combined-spectra analysis relies on the signal adding 
coherently while the noise adds randomly. Thus it is most 
effective if the signal remains in the same bin from spec- 
trum to spectrum. Because the HR channel is affected 
by the Doppler shift of an axion signal due to the time 
varying velocity of the detector, this constraint imposes a 
lower limit on the bin width of this analysis. In a worst- 
case scenario, the combined daily signal modulation from 
both the orbital and rotational motions of the Earth is 
at most a few hertz [25l |26{ . A bin width larger than 
twice the modulation amplitude is required to minimize 
the probability of the signal drifting into a neighboring 
bin between scans. Therefore, each neighboring set of 
256 bins in the 42 mHz spectrum was co-added, yielding 
a spectrum with a bin width of 10.8 Hz. The resulting 
spectrum is sufficiently averaged so as to possess Gaus- 
sian noise and to be mostly insensitive to the signal mod- 
ulation caused by both orbital and rotational terrestrial 
motion on the time scale of a few days. 

The modulation- insensitive spectra, shown in figure |31 
were corrected for systematic effects, the most notice- 
able being the shape given to the power spectrum by the 
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FIG. 3. An example of a modulation-insensitive power spec- 
trum after the crystal filter shape and amplifier chain effects 
have been divided out. The insert shows the raw power spec- 
trum with a passband shape imparted by the crystal filter. 



crystal filter. A reference spectrum representative of the 
crystal filter shape was created by averaging several days 
worth of data. This reference spectrum was divided out 
before each spectrum was cropped to the 30-kHz band- 
width of the crystal filter. Any lingering broad spectral 
structure was due to frequency dependent interactions 
within the amplifier chain. A sixth-order polynomial was 
fit to and divided out of each spectrum to remove this 
structure. 

Frequencies with excessive power (i.e., bins containing 
power which is comparable to the expected power deposi- 
tion from KSVZ axions) and/or excessive noise (i.e., SNR 
below 5) were rescanned within a day or two. These res- 
cans were averaged with the previous data set to reduce 
noise caused by statistical fluctuations. For frequencies 
having too much power to exclude a KSVZ axion signal 
at 90% confidence, rescans totaling an additional 10 to 
20 minutes of exposure time were taken to determine if 
the signal warranted further consideration as an axion 
candidate. All rescans extended above and below the 
frequency in question by several kHz. Because an actual 
axion signal must be persistent, the disappearance of a 
signal would indicate that its origin is either from sta- 
tistical fluctuations or from an external transient source, 
rather than from an axion flow. In the frequency range 
of 812.0 MHz to 892.8 MHz (3.3 ^eV - 3.69 fxeV) we res- 
canned several frequencies with excess power, though no 
signals persisted through the rescan process. 

Because no signals remained after performing rescans, 
the combined SNR and accumulated power were used to 
set limits on the local density of non-virialized axionic 
dark matter. Both KSVZ and DFSZ models are con- 
sidered. In the 800 MHz range, a peak in the power 
spectrum from an axion flow with velocity relative to the 
detector of /3 ~ 10~ 3 and dispersion of A/3 <3x 10~ 6 
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FIG. 4. Density limits for the HR channel (A/3 < 3 x 10" 6 ) 
at 90% confidence from 812 MHz to 892.8 MHz. The scale 
for limits on KSVZ and DFSZ axions are shown on the left 
and right axes respectively. Density limits for the currently 
published MR channel data (A/3 < 
for comparison. 



2 x 10" 4 )[24|] are shown 



would have a width of < 5 Hz. The limits presented 
here are for axion flows having this dispersion or lower. 
At 90% confidence, we limit the local density of non- 
virialized axions to exclude p > 0.2 GeV/cm 3 for KSVZ 
models and p > 1.4 GeV/cm 3 for the DFSZ model. Tak- 
ing 0.52 GeV/cm 3 as an estimate of the average local 
density [27l. [28j . non-virialized axions are excluded for 
KSVZ models over a mass range of m a = 3.3/ieV to 
3.69 ^eV. Consider also the caustic ring halo model, pre- 
sented in detail in Ref. [9j, which predicts that most of 



the dark matter in the local neighborhood is in discrete 
flows. At the Earth's location in the Milky Way, one such 
flow is predicted to have a density as large as p ~ 0.84 
GeV/cm 3 which is well over our stated exclusion limit 
for KSVZ model axions. The HR density limits for non- 
viralized KSVZ and DFSZ axions are shown in Fig. |4] 
with the MR limit (calculated from [24j) overlayed for 
comparison. 

The existing time-series data will permit a future anal- 
ysis of this same data set that will be capable of reaching 
finer resolutions while still accounting for terrestrial sig- 
nal modulation. The analysis will return to a single initial 
scan analysis, where candidate signals found in one scan 
will be compared against any rescans having frequency 
coverage corresponding to where the signal could have 
moved to as a result of Doppler shifts. The time be- 
tween scans will determine the maximum expected sig- 
nal modulation for these candidates. To warrant fur- 
ther consideration, a candidate signal must be present in 
rescans within a window set by the maximum Doppler 
modulation about the original candidate frequency. In 
addition, taking the FFT of smaller intervals of the time 
series data enables a range of dispersions to be tested. 
Further, it can be seen in Eq. [4] that the SNR scales in- 
versely as the square root of the signal bandwith. Thus a 
finer resolution analysis would yield a higher sensitivity 
to g aT y. We expect approximately an order of magnitude 
improvement in sensitivity for the HR channel using such 
an analysis at a resolution of 100 mHz. 
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